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Abstract 

We studied the quantum dynamics of ferromagnetic domain walls (topological kink-type solitons) 
in one dimensional ferromagnetic spin chains. We show that the tunneling probability does not 
depend on the number of spins in a domain wall; thus, this probability can be large even for a 
domain wall containing a large number of spins. We also predict that there is a strong interplay 
between the tunneling of a wall from one lattice site to another (tunneling of the kink coordinate) 
and the tunneling of the kink topological charge (so-called chirality). Both of these elementary 
processes are suppressed for kinks in one-dimensional ferromagnets with half-integer spin. The 
dispersion law (i.e., the domain wall energy versus momentum) is essentially different for chains 
with either integer or half-integer spins. The predicted quantum effects could be observed for 
mesoscopic magnetic structures, e.g., chains of magnetic clusters, large-spin molecules, or nanosize 
magnetic dots. 

PACS numbers: 75.10.Jm, 03.75.Kk 
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I. INTRODUCTION 



Domain walls play an important role in the physics of magnets. For macroscopic bulk 
magnetic samples, domain walls,being extended classical objects, determine the demagne- 
tizing processes, see, e.g., Refs. 1I2J3I. For one-dimensional magnets, domain walls (kink- type 
solitons) play a different role: they are nonlinear excitations responsible for the destruction 
of long-range order. Thus, domain walls should be taken into account together with the lin- 
ear excitations (magnons).- 1 ^ Quantum properties are inherent to kinks in one-dimensional 
magnets (spin chains) with small spin values, like S'=l/2orS=l and high anisotropy.- 

Classical solitons in one-dimensional Heisenberg ferromagnets have been investigated in 
detail. For continuum media, their dynamical properties are determined by the Landau- 
Lifshitz equation for the magnetization vector m(x,t), where m 2 = 1, see e.g., Refs. [j, 
For such systems, kink-type solitons can be treated as classical particle-like objects. 
However, kinks are extended objects, and for spin chains with low anisotropy K <C J (J 
is the exchange integral and K the anisotropy constant), a kink involves a large number 



of spins -/Vjjink ~ SyJ J/K ^> 1. For this reason, domain walls for low- anisotropy magnets 
should be formally considered as a mesoscopic, rather than a microscopic, object. Thus, it 
is not obvious whether or not the effects of quantum coherence are essential for domain wall 
dynamics in mesoscopic ferromagnets. 

Artificial quasi-lD mesoscopic materials (including chains of small magnetic elements, 
such as small magnetic particles of nanometer size (magnetic dots), patterned magnetic 
films, magnetic clusters and high-spin molecules) are of great importance,-&*^>^^^ 
and are promising elements for computers.- These materials often manifest unique physical 
properties that are absent in bulk samples, for instance, macroscopic quantum coherence 
and quantum tunneling, see e.g., Refs. Qfl These quantum properties of small magnetic 
systems could be potentially useful for quantum computing.— 

Quantum coherence can occur when states with the same energy are separated by a small 
potential energy barrier. Kink-type solitons can demonstrate a rich variety of different quan- 
tum effects. Indeed solitons are particle-like objects and their quantum dynamics include 
the tunneling of the kink coordinates through a potential barrier, separating equivalent po- 
sitions in ID chains. These tunneling effects are coherent, and, as for electrons in crystals, 
they lead to the formation of band spectra. 
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Another type of tunneling involves the domain wall chirality. Namely, a domain wall is 
characterized by the deviation of the magnetization from the easy axis. The corresponding 
value of the total spin Sj_ is quite large, i.e., |Sj_| is of the order of SN kink ^> 1. For 
biaxial magnets, the domain wall state has two-fold degeneracy along the Sj_ direction. 
From a mathematical point of view, the sign of the quantity corresponds to the value 
of topological charge or chirality x = ±1. Thus, the kink structure is doubly-degenerate 
over the sign of this topological charge, implying the possibility of a quantum coherent 
superposition of two states with different chirality. 

The tunneling of topological charges has previously been discussed for different topological 
solitons (kinks, vortices, and disclinations) in antiferromagnets; for a review see.— It is 
worth noting here that the static distribution of the corresponding order parameters, the 
normalized magnetization vector m for ferromagnets (m 2 = 1) and sublattice magnetization 
vector 1 (l 2 = 1) for antiferromagnets, are all identical. For antiferromagnetic spin chains, 
the rate of the chirality tunneling process appears to be unexpectedly high, because the 
tunneling exponent is of the order of the atomic spin S, being independent on the number 
of spins iVkink 1 within the kink.— Thus, quantum effects could be essential not only for 
literally ID objects like spin chains, but also for mesoscopic antiferromagnetic samples, like 
thin antiferromagnetic wires.— 

The dynamic properties of ferromagnets are described by the Landau-Lifshitz equation 
for the vector m with no inertial term. This is in strong contrast with the inertial dynam- 



ics of antiferromagnets, described 
magnetization vector 1, e.g., Refs. 



by the so-called sigma-model equation for the sublattice 
3j. It might look a bit paradoxical that the quantum 
properties of ferromagnets are much more complicated than for both quantum antiferromag- 
nets and quantum Josephson junctions^* 2 ^!. The reason is because within the sigma-model 
approach the dynamics of the vector 1 is similar to that of the usual inertial dynamics of a 
particle (strictly speaking, the dynamics of a particle along the surface of the sphere l 2 = 1); 
whereas the Landau-Lifshitz Lagrangian contains a complicated Dirac-monopole term with 



23. This circumstance leads 



non-trivial topological properties (Berry phase), see, e.g., Ref. 
to a number of subtle and intriguing effects, e.g., the suppression of tunneling transitions 
due to the interference of the instanton trajectories.—* 2 ^ 

An example of strikingly different quantum dynamics of ferromagnets and antiferromag- 
nets is the tunneling chirality. For kinks in antiferromagnets, chirality tunneling is not 
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correlated with the translational motion of kinks.— ^ In contrast, some results 
the literature imply that for ferromagnets the situation can be different. Ref. 



known m 



27 noted 



that unmovable kinks with different values of the chirality must have different values of the 



momentum. Ref. 



28 



pointed out that the chirality tunneling rate grows with the intensity 
of the kink spatial pinning. Ref. |29| showed that for a free (i.e., with no pinning) domain 
wall in a ferromagnet, the momentum of which is conserved, the tunneling of the chirality 
is prohibited. 

In this article we develop a consistent quantum theory of domain walls in one-dimensional 
ferromagnets, with a complete description of all possible coherent quantum effects. We show 
that the tunneling of the chirality can be described as a tunneling in momentum space. This 
is closely connected with the tunneling of the kink coordinates; indeed, all these tunneling 
processes are naturally described within phase plane (X, P) . 

The article is organized as follows. In the next section [Til we introduce the Hamilton 
variables: kink coordinate X and kink momentum P. In this section [TTJ we show how to 
consistently define the chirality of a kink via the value of kink momentum. A Hamiltonian 
approach, valid for describing both quantum effects, tunneling of kink coordinates and kink 
chirality, will be developed in the same Section [TTJ Then, in Section IHIt specific tunneling 
effects will be analyzed based on this approach. Kink dispersion relations will also be derived 
there. 



II. HAMILTON DESCRIPTION OF KINK DYNAMICS 

In order to describe a one-dimensional system of mesoscopic magnetic particles allowing 
kink dynamics, we assume that each particle has an internal magnetic anisotropy, with the 
chosen axes to be parallel for all particles in a system. The geometry of the problem is 
shown in Fig. [TJ We will also consider isotropic nearest neighbor interactions. Moreover 
we assume that any internal degrees of freedom of the particles can be neglected, and each 
particle can be treated as a single magnetic moment (spin). 

Note that kinks in a ferromagnet with pure axial symmetry cannot move because the 
projection of the total spin on the easy axis S^ ot ^ should change while the kink is moving. 
However, the Hamiltonian is invariant under rotations around the easy-axis and it commutes 
with S^ ot \ prohibiting such dynamics. In order to allow kink dynamics, we can consider the 
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FIG. 1: The geometry of the problem; a) anisotropy axis for a single magnetic particle; b) a 
schematic representation of a chain of particles. Each particle is represented by its spin. 

two-fold magnetic anisotropy. This is in contrast with antiferromagnets, where kinks can 
move even for purely uniaxial anisotropy, with an easy axis of symmetry. 



A. Model 



The Hamiltonian of a chain-like system of magnetic particles can be written in the same 
form as for a discrete ferromagnetic chain (one-dimensional lattice with atomic spacing a) 
with a spin operator S ra located in each lattice site n, 

K = - J Sn ■ S n+1 + [Ki S tn + K 2 SIS (1) 

n n 

Here the first term describes the isotropic exchange interaction of spins, and the second 
sum corresponds to the two- fold magnetic anisotropy. We have chosen K 2 > K\ > 0, so 
that the orthogonal axes 3, 1 and 2 are the easy axis, the medium axis and the hard axis, 
respectively. We consider each spin operator to be a classical vector with constant modulus 
S = Sm and unit vector m (i.e., this vector points in the unit sphere m 2 = 1). The dynamics 
of the variables m n is governed by the discrete version of the Landau-Lifshitz equation, 

«9m n / dW \ 

where W = W(m n ) is the energy of the considered ferromagnetic chain. The dynamics of 
the variable m n , for a given point n in a lattice, is determined by the Lagrangian, 

£ = -hSY- ^ r • (m n x ^) - W, (3) 

^ 1 + (n ■ m n ) V 9t J 



where n is an arbitrary unit vector. For the continuum approximation, the set of variables 
m„ determined on the lattice sites n should be replaced by a smooth function of the con- 
tinuous coordinate x: m n — > m(x). In this continuum approximation, the Lagrangian takes 
the form 

J \ dt J ' a l + (n-m)' v ' 

where W is an energy functional written by expanding over gradients of m; see, for example, 



Ref. 



24 The vector A has the form of a vector-potential of a Dirac monopole field, i.e., the 



curl of A over the variable m is proportional to m: 

curl m A = HSm/a; 

see, for example, Ref. |23 



The vector potential has a singularity (Dirac string) for (m • n) = 
— 1, i.e., on a half-line in m— space. It is important to note that the vector potential A is 
accurate within some gauge transformations, which includes changing the direction of the 
Dirac string, but the Landau-Lifshitz equations containing rot m A are invariant with respect 
to gauge transformations.— It is worth mentioning here that for antiferromagnets within the 



sigma-model approach, the dynamical part of the Lagrangian has a standard inertial term 

£dyn,AFM OC 



0P 2 



m 

in contrast to A ■ (dm/dt) for ferromagnets. 

It is convenient to represent the unit vector field m by two independent angular variables 
9 and <p, 

mi = sin#cos</9, rri2 = sin 9 sin <p , m^ = cos9. (5) 

In terms of these variables, the Landau-Lifshitz equation takes the form 

Sh . „ 89 5W 

— sm9— = — , 6 
a dt dip 

Sh . n dtp 5W 

— sm6 -£ = - —, 7 
a dt 69 

where W = W{9, <p} is the energy functional written it terms of the field variables 9 and tp. 
The Hamiltonian (TjQ) can be rewritten as 



W{9,^} 



dx Ja 2 



a 2 



09 \ . 2 Jd^ 
+ sin 9 ' 



dx J \dx 



sin 2 9 [K t + (K 2 - K t ) sin 2 <p] \ , (8) 
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where the non-zero difference (K% — K\) determines the anisotropy in the basal plane. 

The Landau-Lifshitz equations have an obvious integral of motion, the energy W. In the 
continuum approximation, translational invariance leads to the conservation of the linear 
momentum of the magnetization field P (momentum, for short). The expression for the 
momentum is determined^ by the dynamical part of the Lagrangian 



As noted above, the vector-potential A is known up to a gauge transformation, and the 
momentum P depends on the gauge used for A. This is one of the problems for describing 
the dynamics of solitons in ferromagnets. This problem does not exist for antiferromagnets, 
where the momentum is proportional to the integral J dx(d\/dx)(dl/dt). However, it turns 
out that the relative momentum for any pair of kinks in a ferromagnet can be uniquely 
determined. Thus, the kink momentum is accurate within the position of the origin in 
P— space; that is, the choice of a kink assigned with the value of P = 0. 

B. Topological analysis of domain wall structure 

For a kink in a ferromagnetic chain, the values of the on-site variables Ss jn have opposite 
values in front and behind of the kink: £3^ — > ±S at n — > ±00. In other words, a kink 
can be seen as a path connecting the poles of the sphere m 2 = 1, corresponding to the 
two easy directions of magnetization space (see Fig. |2j). For definiteness, we assume that 
1713 = 1 at x — > —00 and 7713 = — 1 at x — > +00 (see Fig. EJ). For ferromagnets with non-zero 
anisotropy in the basal plane, dW/dip 7^ 0, a domain wall can move with some velocity v 
smaller than the limit value v c . Within the continuum approximation, such moving domain 
walls are described by a simple traveling-wave one-soliton solution of the Landau-Lifshitz 
equation (jSJ) of the form 9 = #(£), (p = <p(£), with £ = x — vt. However, it is hard to find 
an analytical solution of a set of two second-order equations of this type, and we will start 
with a qualitative analysis. 

Due to symmetry, there are two types of stationary domain walls, with 9 = 9(x) and 
<p = const. For these domain walls, the vector m turns either within the easy plane (3, 1), 
for the first type of walls, or within the hard plane (3, 2), for the second one. The trajectories 
describing these domain walls are denoted by symbols B+, B— and Ni, N2, respectively, on 




(9) 
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FIG. 2: (a) Trajectories on the sphere m 2 = 1 (the direction of motion is indicated by the arrows) 
corresponding to different kinks; the labels 1, 2, and 3 indicate the anisotropy axes of the magnet. 
Two of the most favorable kinks with chirality x = ±1 ar e shown by the symbols B+ and B— ; 
unfavorable kinks with indefinite chirality are labeled by N\ and N 2 - The crossings between the 
trajectories and the axes 1 or 2 are shown by four gray ovals in (a). The trajectories N\ and N 2 
divide the sphere in two domains, associated with kinks of different chirality, as discussed in the 
text, (b) Spin distribution for favorable kinks with chirality % = ±1. 

the Fig. [2j Their energies are E\ for B± and E 2 > E\ for iV lj2 . Other domain walls having 
v 7^ and energies E 1 < E{y) < E 2 , are described by paths, located in between these chosen 
trajectories on the sphere. 

The kink momentum is the total momentum of the magnetization field, calculated along 
the corresponding solution of the Landau-Lifshitz equation.-^ For domain walls, it can be 
written in the form of an integral around the contour of the sphere Fig. [2} depicting a 
kink, P = J A(m) dm, there A is the vector-potential of the Dirac monopole field (j3J). A 
difference of the momentum values for two different kinks can be described as an integral 
along a closed contour. It can be written through a surface integral of the type of J dS rot m A 
and it is equal to hS/a, multiplied by the area, on the sphere m 2 = 1, inside two trajectories, 
corresponding to these pairs of kinks.— It is clear that for a biaxial ferromagnet there are 
pairs of diametrically opposite trajectories (e.g., the trajectories B+ and B— in Fig. [2]) 
corresponding to energetically equivalent but physically different kinks. For these pairs of 
kinks, the closed path borders half of a sphere, with an area of 2-ir, and the momentum 
difference equals 2nhS/a. Hence we can readily obtain the periodic dependence of the kink 



S 



energy on its momentum with the period Po, 

„ 2ttHS 



o 



(10) 



a 

For those kinks which are close to the most favorable kink B+ (B-), the value of the 
momentum P (AP = P — P ) is small; thus, the parabolic approximation can be used and 
the energy can be written as E = P 2 /2M (or E = (AP) 2 /2M), where M is the effective 
mass of the kink. For the model (JSj) , the effective mass M takes the value of the well known 
Doring effective mass, obtained as early as the 1930's; see, e.g., Refs.|2|j6|. This effective mass 
turns to infinity when {Ki — Ki) — > 0. This is another indication that in a pure uniaxial 
model of a ferromagnet (K 2 = Ki), domain wall motion is impossible. However there is no 
contradiction between the finite value of P oc (p and the condition v = 0: if P = Mv , then 
the momentum can be finite when M — > 00 and v — > 0. 

Coherent tunneling assumes the presence of at least two different states having the same 
energy; for instance, the two states of kinks in a biaxial ferromagnet with different values of 
the topological charge. A topological classification^^ of kinks can be done in the same 
way for both ferromagnets and antiferro magnets. First, the difference of values of the 
magnetization vector (or l x , for antiferromagnets) on the right and on the left sides of 
the kink determine the topological charge 7To of the kink. Changing this topological charge 
requires overcoming the potential barrier, proportional to the system size (formally, infinite 
barrier) that cannot be realized by tunneling. Second, a topological charge of the type 7Ti 
is determined by mapping the coordinate space of the spin chain (the line — 00 < x < 00) 
onto the circle {m 2 + m\ = l,m 2 = 0}, situated in the easy-plane of the ferromagnet. The 
appearance of two topological charges of different levels can be formally described using the 
relative homotopy group, as discussed in Ref. 
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The 7Ti charge is described by the integral 



X = - e 2 m x — — dx. 



TV 



dx 

In other words, the chirality x — ±1 determines the sense of rotation (clockwise or coun- 
terclockwise) of m along the chain. This standard definition of the chirality used in 



Refs. 
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I l29l is only valid for kinks located on the unit sphere near the most energet- 
ically favorable configurations B+ and B-, as shown in Fig. [2] that is, for kinks having 
small velocities v. For this case, the effective mass approximation is valid, and these two 



kinks are well separated. According to the topological analysis, these kinks have different 
values of the chirality \ = ±1; and within a self-consistent Hamiltonian approach, they 
correspond to different values of the momentum: P = and P = Pq. However, for kinks 
moving with a non-small velocity and with arbitrary values of the momentum, the above 
definition of the chirality should be modified. 

For treating the whole order parameter space (i.e., the sphere m 2 = 1), kinks with x — ±1 
can be transferred to each other through energetically unfavorable kinks of type N\ or N2, 
schematically shown in Fig. [2l Here the barrier is finite (and equal to E 2 — Ei\ which is large 
when K2 ^> K{) and the process of kink chirality x tunneling is possible. The concept of 
chirality, as a discrete number x — =tl? is naturally connected with the presence of a discrete 
degeneracy in the dependence of the kink energy on its momentum. The discrete parameter 
X = ±1 determines one of two different, but energetically equivalent, kink states existing 
in a biaxial ferromagnet. The values x = 1 an d X = ~ 1 can De naturally attributed to kink 
states with trajectories in two equivalent semi-spheres, m\ > and mi < 0, respectively. 
The chirality value is not determined for the unfavorable static kinks only (Neel walls) 
for which the trajectories N\ and N2 pass through the hard axis. In this sense, chirality 
tunneling can be seen as a tunneling effect in momentum space with a non-small (of the 
order Pq = 2irhS/a) change of the kink momentum. 



C. Moving domain wall structure 

To confirm the general features for moving domain walls discussed above, we will discuss 
an exact solution of ([6]) known for the model of biaxial ferromagnets with the energy (jSj). 
It is easy to find the structure of domain walls with zero velocity. There are two types 
of domain walls having thickness = a a/ J/2K\^ and energies E\p = 2Sy/2JKi^ with 
Ei < i?2- These are the one-dimensional analogs of the usual Bloch and Neel domain walls; 



• 3:6 



see Refs. 

The structure of a domain wall moving with a non-small velocity within the model (JSJ) 
was obtained by Walker at the end of the 1950s, see Refs. For this solution, the 

function 9 = £ = x — vt, and the value of (p — (po — const is independent on £. The 
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value of ipo is determined by the domain wall velocity v as follows 



vh e sin ip cosy? 

Q —n . =» (11) 



a\J JK + e sin 2 ' 

here and below we use the notation e = {Ki — K\)jK\ to shorted the expressions. The 
relation (ITTj) governs, in particular, the maximal possible value of a domain wall velocity, 
the so-called Walker velocity v w 

v w = (a/h)^fjK(VTT^ - 1). (12) 

It is interesting to note that vw is smaller than the minimal phase velocity of spin waves, 



v ph = (a/H)VJK(VTT^ + 1). 

The values of vw and t> p h coincide only in the limit e — » oo, when the Landau-Lifshitz 
equation (jSJ) can be mapped onto the sine-Gordon equation. However, if e is small, then 

v w - ^ph/4 < w ph - 

The value of vanishes when e —>■ 0; that is, as was noted above, the domain wall cannot 
move at all for pure uniaxial ferromagnets. 

The Walker's solution can be presented in the explicit analytical form 

r £ i r r ^ n- 1 

cos 6 = q ■ tanh 



x {v) 



sin 9 = a < cosh 



(13) 



where g = ±l,a = ±l, 



xn(v) = ail ; r — - (14) 

U 2 J ff 1 (l + £sinV) v ; 



and ip = <p(v) is determined by Eq. (TTTT) . The two topological charges introduced above 
naturally appear here as the quantities q and a; q = ±1 determines the 7r — topological 
charge of the kink, and a = ±1 governs the spin direction in the kink center and it is 
naturally connected with the 7Ti —topological charge, the chirality x- 
A straightforward calculation of the kink energy leads to the formula 



E — Exyl + e sin 2 <p, E x = 2S 2 y / 2JK 1 , (15) 

where <p = <f(v) is determined by Eq. (Ilip . The E(v) dependence consists of two branches. 
In three-dimensional ferromagnets, where two-dimensional plane domain walls are present, 
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the upper branch of the domain walls is unstable. However this instability is developed via 
non-uniform perturbations in the domain wall plane, and such fluctuations do not exist for 
domain walls in one-dimensional magnets. Below we will show that for the more natural do- 
main wall energy representation (namely, as a function of its momentum) the E(P) function 
is single- valued. 

For the Walker solution, the value of (p is ^—independent,- 1 ^ and the paths on the sphere 
present at Fig. [2] are the halves of the big circles passing through the poles of the sphere. 
Therefore, the domain wall momentum can be written as 

P = Po~, (16) 

71 

and the explicit form of the dependence of the domain wall energy on its momentum can be 
rewritten as 

f 2T . 2 nP 2T K 2 - g. 



where 



Ei = 2S 2 y/2JK~i 



is the minimum energy of the domain wall introduced above. It is worth noting that the 
equation (|T7j) reproduces the dependence v(<p) ffTTj) within the Hamilton dynamics, 

dX = dH(X, P) 
~dT =V ~ dP ' 

This leads to unusual dynamical features, like the oscillatory motion of a domain wall as 
a response to a dc driving force (e.g., dc magnetic field parallel to the easy axis). These 
peculiarities are well-known for the exact Walker solution and have been experimentally 
established for moving domain walls in magnetic bubble materials; see, e.g. Ref. [J Here we 
were able to write down the explicit form of the function E(P), but the periodic dependence 
E{P) with the same value of Pq is present for any continuum model of a ferromagnet with 
biaxial anisotropy. Such periodic dependence E(P) is also valid for biaxial discrete models; 



for details, see Ref. 
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D. Kink coordinate and lattice pinning 

As we will show, the quantum properties of kinks can be described within a semiclassi- 
cal analysis of the Hamilton dynamics of collective variables: the kink coordinate X and 
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conjugated momentum P. This dynamics is determined by the characteristic Hamilton func- 
tion H(P,X). In a continuous approximation, the definition of kink coordinate is obvious. 
However, the Hamilton relation 

dP dH(P, X) 

~dt ~ OX 

shows that the kink momentum is conserved for any model with the Hamiltonian independent 
on X. Therefore any processes of tunneling in momentum space (in particular, the tunneling 
of the domain wall chirality) can only occur if the Hamilton function H depends on the 
domain wall position X. For our model, such dependence can only be caused by a lattice 
pinning of the kink. Thus, for a consistent description of quantum tunneling, lattice pinning 
must be considered. 

A first step in this direction is to define the domain wall coordinate X treated as a 
collective variable and conjugated to the kink momentum P. The kink coordinate X in the 
discrete model can be naturally determined through the spin operators, 

n 

, where corresponds to a certain "reference" kink=, which coordinate is chosen as 
X = 0. The total spin projection onto the axis 3 is conserved, Sl ot = $3 = const; thus, 

n 

for uniaxial ferromagnets with K2 = K\, dX/dt = and a kink dynamics is impossible. 
Another consequence appears when taking into account the Hamilton relation dX/dt = 
dH(P,X)/dP. For the purely uniaxial case, H(P,X) does not depend on P. All of these 
general considerations here are characteristic of the exact Walker solution. 

A spin configuration corresponding to a kink with a specified value of the coordinate 
X can be obtained by minimizing the Hamiltonian with respect to the variables 9 n and 
(p n , for a fixed value of the total spin S^ ot \ To do this, we use a procedure proposed and 
numerically realized for the analysis of different dynamical solitons,— which are described 
by a conditional minimum of a discrete spin Hamiltonian. Using this method one can easily 
determine the structure of the kink and obtain the dependence of the kink energy on its 
coordinate X; this for finite spin chains described by any classical spin Hamiltonian. 

Now, we consider model (pQ) with a purely single-ion anisotropy, for a finite chain of size 
N c , with boundary conditions cos8 n = 1 and cos# n = —1 at different ends of the chain. The 
size of the chain N c is chosen to be much larger than the width of the kink. In particular, 
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FIG. 3: Shape of the domain wall pinning potential U(X), normalized by its maximal value U(a/2), 
for two values of the anisotropy constant. Symbols denote the numerical data. The fit for the model 
dependence (|19p is shown by the full line. 

for a reasonable anisotropy K > 0.2J, the kinks occur as well-localized excitations. As a 
result, the kink energy is independent of N c for N c > 30. For extremely high values of the 
anisotropy, K > K c , with K c = 0.667J, the domain wall becomes purely collinear,— with 
all spins up or down, 5*3 = ±S. For such collinear states, the continuum description of the 
domain wall dynamics and its topological analysis are obviously incorrect, and we should 
restrict our consideration to moderate values of the anisotropy: 0.2J < K < 0.65J < K c . 

Considerable influence of the lattice pinning appear when K > 0.25J. A more favorable 
position of the kink is between neighboring spins. Thus, the values Ss jtl at two neighboring 
spins are equal in magnitude and opposite in sign: S3 = ±5(0) < 1. Choosing the value 
of X = for one of such states, we can determine the pinning potential having equivalent 
minima at the points X = an, where n is an integer. The states with kink on a lattice site 
with X = a(2n + l)/2 correspond to maxima of the pinning potential, as shown in Fig. [31 
In general, it can be concluded that, for moderate values of the anisotropy K < K c , the 
pinning potential is not large compared to the "static" energy of the kink. The dependence 
U(X) is fairly well described by the simple harmonic relation 



On the other hand, higher Fourier components are also present in the dependence U(X), 




(18) 
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FIG. 4: Coefficients Uq and U\ (in units of JS 2 ) in Eq. (fl~9|) for some values of anisotropy constant 



the contribution of U\ can be noticeable. The values of Uq and U\ as functions of the reduced 
anisotropy Kj J are presented in Fig. |H In general, for higher values of the anisotropy K 
one can see a broadening of the curve U(X) near its maximum and, correspondingly, a 
narrowing of the curve near the minima. Note that when K > K c , the function U(X) has 
a cusp at X = 0. 

To conclude this section, we derive the periodic dependence of the domain wall Hamilto- 
nian on both collective variables X and P. Namely, to describe the quantum dynamics, we 
can use the Hamiltonian H(P,X) = T(P) + U(X), where both functions T(P) and U(X) 
are periodic: U(X) = U(X + a) and T(P) = T(P + P ). 

III. QUANTUM TUNNELING EFFECTS IN KINK DYNAMICS 

To describe the quantum dynamics of domain walls, we can proceed with the Hamiltonian 
H(P,X) = T(P) + U(X), where the "kinetic energy" T(P) is described by Eq. (02} and 
the "potential energy" is caused by a periodic lattice pinning potential U(X) = U(X + a); 
see Eqs. ({IB}) or (fl9l) . The most crucial point is the presence of the double periodicity of 
H(P,X), with respect to both X and P. For simplicity, below we use the simplest model 



especially for higher anisotropy. For example, for a more general form 




(19) 
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holding this property 

if = Tosin 2 ^+f/ sin 2 ^, (20) 

where only the lowest harmonics are considered. Here the energy of the kink is taken from 
its minimum classical value, i.e., the energy of a "unmoving" kink with P = or P = P , 
situated at the minimum of a pinning potential X = 0. 

For different ferromagnetic chains, any ratio of parameters T and Uq is in principle 
possible. For example, the value of T is proportional to the difference (K 2 — Ki) and it 
vanishes in the purely uniaxial case for any value of K 2 — K\ ^ 0, while the amplitude 
of the pinning potential Uq is almost zero for K\ < 0.15 J. It is natural to start with the 
Bohr-Sommerfeld quantization for the domain wall motion, which is based on the analysis 
of the classical dynamics. This can be done in the same way as for the transverse-field Ising 
model (see Ref. 1221 ) and we do not discuss its details here. The most important feature of 
a Hamiltonian of the type (|2Q|) is the presence of a lower and upper bound for the energy. 
Hence, two types of finite motion appear. A first type corresponds to oscillations of the 
domain wall with either P <C P or (P — P ) <C P near the minimum of the potential, 
with the energy near the minimum of the Hamiltonian, E <C To, Uq. A second type of finite 
motion corresponds to oscillations near the potential maximum, with the momentum near 
the values of P = ±Pq/2, the values of energy are E < T + Uq. 

For intermediate values of the energy, the motion is infinite. For a small pinning potential, 
this motion is standard, with infinite grow of the kink coordinate. The case Uq > T is 
less standard; it corresponds to an infinite-growing momentum with finite oscillations of 
the coordinate near certain positions, which do not coincide with extrema of the pinning 
potential. The late case is nothing but Bloch oscillations in the pinning potential U(X). An 
exception is the chosen point To = Uq, for which the classical motion is finite for all values 
of the energy. 

Both types of infinite motion describe the classical over-barrier dynamics of the domain 
walls. Using a quantum-mechanical language, such states of the nearly-free particle can be 
well described by perturbation theory over Uq or T , for the cases Uq < T or Uq > T , 
respectively. This analysis will be done in the next subsection IHI Al 

Within the Bohr-Sommerfeld quantization condition, all the states corresponding to finite 
motion (oscillations) of the domain wall near any extrema of the Hamiltonian produce a 
discrete spectrum of energy levels E n with a level separation of the order of \/UqTq/ S. Both 
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Uq and T are proportional to S 2 , therefore for the semiclassical situation of high spins 
S > 1, the separation of values of E n oc S and it can be smaller than Uq or To. These 
states, with energy E n , are well localized. For them, the probability of tunneling is small, 
and for its estimate the semiclassical approximation is adequate. Such analysis will be done 
in subsection IIII Bl 

A. Pert ur bat ive analysis 

For extremely large or small values of U /T (namely, for U /T < 1/S 2 or T /U < l/S 2 ) 
the value of E n+ i — E n ~ fkj n can be of the order of min[T , Uq]. In these cases, all the states 
are delocalized, and our perturbation theory gives the full description of the domain wall 
spectrum. Out of these strong inequalities, a perturbative analysis can only be applied to 
domain wall states with intermediate values of the energy, which correspond to the classical 
infinite motion discussed above. 

Let us now start with the case Uq <C To, using perturbation theory with respect to U(X). 
In this case, in zeroth approximation, P = const. To proceed further, let us assume the 
chain to have a large but finite size L = Na, N ^> 1. Then periodic boundary conditions 
give the usual quasi-continuous spectrum of the momentum P = P n = (2tt h / a) (n/N), where 
n is an integer, N/2 < n < N/2, or —Pb/2 < P < Pb/2, where Pb = (2nh/a) is the size 
of the usual (crystalline) Brillouin zone. The vector of states \P) corresponds to a fixed 
value of momentum and, hence, the fixed value of chirality and the uncertainty value of 
the kink coordinate. The quantum spectrum of the problem repeats the dependence of the 
Hamiltonian (I20I) on P. 

We now consider the term U(X) as a perturbation. Its role will generally be the same 
as for the lattice potential U(X) in the standard weak-binding approximation in solid state 
physics. At zeroth order approximation in the coordinate space, the eigen-functions are of 
the form ip^ = exp(iPX), with the energy E(°)(P) = T(P). The influence of the potential 
U(X) with the period a leads to the formation of Bloch states which are a superposition of 
the states ip^ n ' = exp(iPX + inPs), and the momentum transforms to quasi-momentum. 
In the weak-binding approximation, the spectrum can be obtained by a superposition of 
unperturbed dispersion curves E^ = T(P), with argument shifting by nPg, where n is 
an integer number. This spectrum is periodic with the period equal to the size of the first 
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Brillouin zone Pb- The influence of the perturbation is maximal if the values of the functions 
(P + nPs) and (P + ri'Pg), with the different n ^ n', coincide for some value of P. 

In contrast with Bloch electrons with a parabolic dispersion law E^°' el ' = P 2 /2M, for 
kinks in ferromagnets, the unperturbed dispersion law is already described by a periodic 
function. Hence, for the resulting dispersion law E(P), the periodic dependence (with the 
period matched with both characteristic values Pq = 2irSh/a and Pb = 2nh/a) should 
appear. It is also important that kinks for the states \P) and \P + Pq) have the same energy 
and velocity, but differs by the sign of the chirality x — ±1- I n pure classical language, these 
states are described by different magnetization distributions, their images corresponding to 
diametrically opposite paths on the sphere in Fig. [2J 

A simple analysis shows the fundamental difference between the character of the spectrum 
for integer and half-integer values of the atomic spin S. For minimal integer S — 1, the 
periods P and Pb coincide: see Fig. [5j The accounting of the potential U (X) of the form (|20|) 
leads to the overlapping of functions E^°'(P) and E^°\P + Pq). Taking, for definiteness, a P 
situated in the first Brillouin zone, —Pq/2 < P < Pq/2, we can say that these unperturbed 
states have a different chirality x = +1 an d X = ~ 1- The action of the potential leads to 
their hybridization and formation of the states |P±) = (\P) ± \P + P ))/\/2, having energies 

E&\P) = E®\P) ± (P\U{X)\P + P ), 

where 

E®(P)=T(P) + (U(X)), 
(P\U(X)\P + P )) = U /4 

and 

(U(X)) = (P\U(X)\P)) = U /2 

are the off-diagonal matrix element and the mean value of the potential U(X) for (1201) . 
respectively. 

The same expression can be found for any integer spins S = k > 1, i.e. the Fourier com- 
ponent of the potential with P = kPs = Po leads to a full hybridization of the chirality for 
any P. As a result, states of type |+) ± |— ) appear, where |±) correspond to chirality values 
X — ±1- Such states are the quantum superposition of the kinks describing diametrically 
opposite trajectories on the sphere and energy £(±)(P) = £(°)(P) ± U /4, where U is the 
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FIG. 5: Dispersion relation of a kink in a ferromagnet with spin S = 1 subject to a weak pinning 
potential. The solid line and the dashed line correspond to the states (|+) — |— )) and (|+) + |— )) 
, antisymmetric and symmetric over chirality, respectively. Here, and in the Figs. [6HH1 the vertical 
dotted lines show the boundaries of the usual Brillouin zone, — Pb/2 < P < Pb/2- 



S = 2 E(P)IT Q 




v pip. 
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FIG. 6: Same as for Fig.[5]for spin value 5 = 2. The region near the crossing points, schematically 
shown by shadowed ellipses on the main figure, are magnified in the insert upward left. 

corresponding matrix element of the potential U(X). For integer spins S > 1, the value 
Pq = S ■ Pb is a common period. Neglecting the chirality tunneling, one can find the S 
usual energy bands (doubly degenerated over the chirality values) with the size of Pb- The 
chirality tunneling splits any of them into two subbands, corresponding to states (|+) ± |— )) 
with E^(P), and the total number of bands equals to 2S, as shown in Fig. El 

For half-integer spins S = k + 1/2, where k is an integer number, the situation is com- 
pletely different. It is easy to show that none of the Fourier components of the potential 
with uPb = 2nP /(2k + 1) leads to such an overlapping of the non-perturbed spectrum at 
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FIG. 7: Dispersion relation of a kink in a ferromagnet with spin S = 1/2, subject to a weak pinning 
potential. The solid line and dashed line correspond to the states with chirality x = — 1 an d X = L 
respectively. Rectangles denote the points where the value of the chirality is not determined. 




FIG. 8: Same as in Fig. [7] for the spin value S = 3/2. The gray and the light circles schematically 
show the areas in the vicinity of the crossing of unperturbed spectra with and without chirality 
hybridization, respectively. 

any P, which take place for integer S. For S = 1/2 any crossing is absent, as shown in Fig. [7J 



27 is an artifact of 



[a crossings of non-perturbed spectra for spin 5=1/2 reported in Ref. 
the parabolic approximation for E^°\P)]. For higher half-integer spins S > 1/2 the only 
crossings at some fixed values occur when P = P n . Such crossings can appear for branches 
E(°>(P + nPs), corresponding to the kinks with the same or different chirality, as shown in 
Fig. [HJ For this last case, the effects of chirality hybridization can be present very near the 
crossing points, \T(P) — T(P n )\ Uq, as shown in Fig. El 

Let us consider the opposite limiting case Uq ^> T , when the kinetic energy T(P) plays 
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the role of a small perturbation for the unperturbed Hamiltonian Hq = U(X). To zeroth 
order approximation with respect to To, the solution is now X = const. To construct this 
perturbation theory, the momentum representation should be used. In this case, taking 
into account the identity of the states with P and P + 2P , we should apply the condition 
exp[i(P — Po)X] = exp[i(P + P )X], that leads to the discreteness of the values of the kink 
coordinate, X = Xk = ak/2S. The difference between this boundary condition and the 
condition which was used for the case Uq <C T above, leads to essentially different results. 

For the case of strong pinning potential U ^> T , unperturbed states are described by the 
wave function in the momentum representation tp^ = exp(iXkP), with definite coordinate 
X = Xj, and indefinite value of the momentum, which also means an indefinite value of 
the chirality. Let us now consider T(P) as a small perturbation. It is easy to see that nth 
Fourier component of T(P), with exp(2iimP / Pq) , leads to a nonzero matrix element of the 
quantum transition when changing the kink coordinate, X^ — > X^ + 5X n , 5X n = na/S 
(in the simplest case with one harmonic only, Eq. (l20l) . the transition with n > 1 requires 
accounting the nth order in perturbation theory). If U(Xf.) = U(Xk + SX n ), i.e. when the 
value of SX n is multiple to the chain period a, the "resonant" transition should be observed. 

Summarizing, kinks in ferromagnetic chains with either integer or half-integer spins be- 
have differently. The tunneling from a cell to the neighboring one is possible for an integer 



spin, while for the half-integer spin only the transitions with changing X 
ted, i.e. "jumps" across one cell. This feature was mentioned in Ref. 
reasoning. 
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3y 2a are permit- 
using a different 



B. Semiclassical dynamics 

For large values of the spin, S ^> 1, a semiclassical approach provides a much better 
approximation than the perturbative approach used above. For electronic states in a large 
lattice potential, the semiclassical tunneling can be formally described by the tight-binding 
approximation, which leads to the formation of a band spectrum of the kink with narrow 
allowed bands. 

Considering possible tunneling effects for a kink, one immediately encounters the ques- 
tion of under-barrier transition in both coordinate space and momentum space. For domain 
walls with the minimal energy, such transitions include the tunneling between states cor- 
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responding to a two-dimensional set of points in phase space of a system (P, X) , such as 
X ~ 0, ±a, . . . and P ~ 0, ±Po, • • • • Within the semiclassical approximation, these transi- 
tions can be investigated in the framework of the instanton approach; see Refs. 
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38. This 



approach is a version of the Feynman path integral method suited to the description of the 
underbarrier transitions. It involves using the Euclidean space-time, that is transforming 
to the imaginary time, t — > ir (so called Wick rotation). Within this approach, the am- 
plitude of the underbarrier transition from a given quantum state \i) to another one |/) is 
determined by the path integral J DX ■ exp[— where DX denotes integration over 
all possible paths that satisfy the specified boundary conditions. Here the Euclidean action 
A^LY] is described in the form Ae = J LecLt, and Le = P(dX/dr) — H(P,X) is obtained 
by the application of a Wick rotation to the usual mechanical Lagrangian. The instanton 
solution determines the trajectory for which the tunneling amplitude is maximal, i.e.; the 
instanton trajectory minimizing Ae with respect to X{r) and P(t), with the conditions 
|z) at r — > — oo and |/) at r — > +oo. The minimum of the Euclidean action is realized 
on the separatrix solution of the corresponding Euler-Lagrange problem for the Euclidean 
action functional Ae or, equivalently, on the solution of the Hamilton equation with the 
substitution t — > ir. The tunneling splitting of the levels A is determined by the formula 
A = rjhwQ^J A E /hexp(— A E /H), where u is the characteristic frequency, r] ~ 1. 

Let us apply this approach to the mechanical problem of the kink dynamics described by 
the Hamiltonian ( J20l) . It is easy to see that the Wick rotation t — > ir, simultaneously with 
the simple substitution X — > zH, reduces the instanton problem to the Hamilton problem 
for real canonical variables 5 and P and with the real Hamilton function He, 

He = To ■ sin 2 (vrP/P ) - U(S), U(S) = U ■ sinh 2 (vrS/a) (21) 



The Hamilton equations for ( 1211) have an obvious integral of motion He = const; boundary 
conditions give H E = 0. Thus, for an instanton solution, we derive a/^o sin(7rP/P ) = 



±v/U(S). A simple analysis shows that this problem has an instanton solution with 
S(±oo) = (i.e., X(±oo) = 0), while the values of momentum differ: P at r — ► — oo 
and P + Pq at r — > +oo. This instanton solution describes the tunneling of the kink chiral- 
ity 

Analysis of the second type of tunneling (tunneling of the kink coordinate) can be done 
in a similar manner, by using a Wick rotation, substituting P — > ill, and keeping the 
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FIG. 9: Dependence on the parameters and He, determining the value of the Euclidean action 
for the tunneling processes, as a function of Uq/Tq (in logarithmic scale). Two analytic asymptotic 
dependencies are also indicated in the figure with a dashed and a dot-dashed curves. 

coordinate X as a real variable. Then, again, the real Hamilton function He for two real 
variables X, IT appears, 



It is easy to show that this problem has an instanton solution with Il(±oo) = 0, X(r — > 
— oo) — X(t — > oo) = a. This instanton solution describes the tunneling of the kink from 
one lattice site to a neighboring one. 

For the instanton solutions of both types He, He = const = the Euclidean Lagrangian 
reduces to L E = P(dX/dr). Therefore for both cases, (l2T!) and (l22|) . the value of the 
Euclidean action can be represented as simple integrals. For example, for the tunneling of 
the chirality A E = J Q Po H(P) dP, where 



where (Ie(U /T ) is a universal function (see Fig. [H]) that only depends on the ratio U /T . 
We can easily find the asymptotic behavior of a^: the functions a(z) = yfz and a(z) = Inz, 
respectively, where z = U /T Q . These are shown in Fig. 9. 



H E = U - sm 2 (nX/a) - T ■ sinh 2 (7rIT/Po)- 



(22) 



E(P) = (a/ir) arcsinh[V(To/C/o) sin(7rP/P )], 



resulting in 




(23) 
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For the analysis of the tunneling of the coordinate X it is not even necessary to calculate 
the corresponding integral A° E = J a Il(X) dX. The Euclidean action Ae, describing the 
tunneling of the kink coordinate, is obtained from the previous formula for the tunneling of 
the chirality A E in by the replacement U /T — > T /U . 

The quantity Ae, as well as the tunneling splitting A , is a standard measure of the 
intensity of the quantum tunneling processes; but there is a physical difference among the 
tunnel processes considered here. This is because all of the values of the kink coordinate of 
the form X = an (the kink is found at points of potential energy minima) correspond to 
different states of the kink, while all the values of the momentum differing by 2Pq correspond 
to the same state of the kink. Therefore, the processes of tunneling either coordinate or 
momentum lead to different physical consequences. Tunneling of the coordinate is responsible 
for the formation of an energy band in which the number of states with different values of 
the quasi-momentum coincides with the number of spins in the chain. The width of the 
corresponding energy band is given in terms of the tunneling splitting Ao, AEq = 2Aq. 
Tunneling of the momentum leads to lifting off the twofold degeneracy of the states of the 
kink and causes a splitting of the level into two, with AE = 2A 

IV. CONCLUDING REMARKS AND RESULT DISCUSSION 

In conclusion, using both semiclassical instanton and perturbative approaches, we study 
quantum tunneling effects in P— space and X— space for domain walls in ferromagnetic 
chains. We also investigated the quantum dynamics of domain walls (kink-type solitons) 
in spin chains. Explicit results have been obtained for the biaxial model with isotropic 
interaction J, and rhombic anisotropy with two constants K\ and K 2 . The combinations of 
the two spin interaction constants, K\j J and (K 2 — Kx)/K\, define two parameters, Uq and 
To, for the effective Hamiltonian describing the quantum dynamics of a domain wall. 

In both limits (small and large Uq/T ) only one type of transition becomes important, 
but for the case Uq ~ T the probability for both transitions are comparable. In this case 
Uq ~ T , the function a(U /T ) is of order of unity, and for S ~ 1 both transition amplitudes 
are not small; i.e., the kinks in ferromagnetic chains with spin S ~ 1 are essentially quantum 
objects. They are characterized by a quantum dispersion relation (spectrum) of the kink 
E = E(P), with the presence of some discrete variable, chirality \ — ±1. The quantum 
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properties of domain walls in a chain with either integer or half-integer spin are essentially 
different. For a chain with integer spin, there are S main energy bands, each one of them 
split in two subbands, with a total hybridization of the chirality. In contrast, for the case 
of half-integer spin chain we arrive at a pattern of 2S nonoverlapping energy bands, with 
chirality hybridization only at some particular points. 

Having in mind the case of mesoscopic chain-like artificial ferromagnetic structures, we 
discuss the behavior of domain walls for large spins S ^> 1. At a first glance, for such 
systems, the quantum effects should be suppressed by the large spin values. However, as we 
have shown, tunneling effects can occurs, even for values like S ~ 10 2 — 10 3 , for essentially 
different values of the parameters Uq and T . Here the value of the tunneling exponent can 
be acceptable (Ae/K) < 15 - 20 if cle <C 1. We stress the agreement between the probability 



37 



38 



and the probability of 



of tunneling for one mesoscopic magnetic particle, see Refs. 
tunneling processes for the kinks found here. 

For domain wall tunneling, at least one of two quantum tunneling transitions (either 
tunneling of coordinates or tunneling of chirality) are possible. It is useful to introduce 
the following empiric rule. The probability of chirality tunneling, including flipping a large 
number of spins 3> 1 has the same order of magnitude as for coherent spin tunneling 
of a single particle. 
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